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The cos 2(f> asymmetry of Drell-Yan and J/ip production in 
unpolarized pp scattering 
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Abstract. We investigate the cos 2<j> azimuthal asymmetry in Drell-Yan and J/ip production from unpolar- 
ized pp scattering at GSI-HESR energies. The contribution to this asymmetry arising from the leading- twist 
Boer-Mulders function hi (a:,ky), which describes a correlation between the transverse momentum and 
the transverse spin of quarks in an unpolarized hadron, is explicitly evaluated, and predictions for the GSI- 
HESR kinematic regime are presented. We show that the cos 20 asymmetry is quite sizable both on the 
J/ip peak and in the Drell-Yan continuum region. Therefore these processes may offer an experimentally 
viable access to the Boer-Mulders function in the early unpolarized stage of GSI experiments. 
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1 Introduction 

The spin structure of hadrons is under intense investiga- 
tion, both theoretically and experimentally. In hadron- 
hadron scattering one would naively speculate that the 
polarization of at least one of the incoming hadrons, the 
beam or the target, or of both, is needed in order to in- 
vestigate the spin properties of hadrons. However this is 
not the case if one takes into account the intrinsic mo- 
mentum of partons. It is known in fact that there is a 
new kind of leading-twist transverse-momentum depen- 
dent distribution function, the so called Boer-Mulders 
function h^(x, ky) [Tj[2]) which describes the transversity 
of quarks inside an unpolarized hadron. This new function, 
which is chirally odd, manifests itself through the cou- 
pling to another chiral-odd quantity, such as the transver- 
sity distribution hi, the Collins function H± , or another 
Boer-Mulders function [3] . Therefore the spin structure of 
hadrons can be also studied in physical processes without 
beam and/or target polarization. It has been argued by 
Boer 2 that h±(x, kip) can account for the large cos 20 
azimuthal asymmetry observed in the unpolarized pion- 
nucleon Drell-Yan process [HH]. The origin of hi (%, k|,) 
has been addressed in Refs. 0(7] • The first calculation of 
hi(x, k;p), in a quark-scalar diquark model, was reported 
by Goldstein and Gamberg [8]. Soon after this work, Boer, 



Brodsky, and Hwang [5], using a similar model, showed 
that the analyzing power v of the cos 2</> asymmetry could 
be as high as 30% in the unpolarized pp Drell-Yan process, 
and smaller in pp scattering. The Boer-Mulders function 
of the pion was calculated [TOlTTTj in a quark-spectator- 
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antiquark model and shown to reproduce, in combination 
with the corresponding distribution of the proton [12| , the 
cos 2(f> asymmetry in the tt~ N Drell-Yan process measured 
by NA10 Collaboration [4]. The Boer-Mulders function, 
combined with the Collins function, can also produce the 
cos 2(f> asymmetries of semi-inclusive pion leptoproduction 
from unpolarized nucleons jl3J14j. 

Recently, the cos 2(f> azimuthal asymmetry of the pp 
Drell-Yan process has received much attention, as there 
have been proposals to study spin phenomena in polarized 
and unpolarized pp scattering at the High-Energy Storage 
Ring (HESR) of GSI [T5l[ro| . There have been numerical 
simulations [17] as well as model calculations [9lfT8] of the 
cos 2(f> asymmetry in the unpolarized pp Drell-Yan pro- 
cess. The purpose of this paper is to investigate the cos 2(f> 
asymmetry of Drell-Yan and J/ip production in unpolar- 
ized pp scattering, using a model of hi(x,\s.^) [TT1 J T21 J T4] 
that has been adjusted to fit the NA10 tt~N Drell-Yan 
data 4J and shown to be in agreement with the ZEUS [19] 
and EMC [3U] semi-inclusive DIS data. As we will see, in 
the GSI kinematic domain the cos 2(f> asymmetry turns out 
to be rather large. In particular, the advantage of study- 
ing the J/ip production is that, while its cos 2</> asymmetry 
is similar in size to the corresponding asymmetry of the 
Drell-Yan continuum production, the counting rate on the 
J / ip peak is two orders of magnitude higher compared to 
the region above this peak. 
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2 The cos 2< 
production 



asymmetry in Drell-Yan and J/ip 



The Drell-Yan process represents an ideal window on the 
hadron structure, since it probes only parton densities. We 
are interested in the transverse momentum distribution of 
the lepton pairs, hence we have to consider the intrinsic 
transverse momenta of partons inside the hadrons. The 
non-collinear factorization theorem for Drell-Yan process 
has been recently proven by Ji, Ma and Yuan [5T] for 
Qt < Q- 

The angular differential cross section for the unpolar- 
ized Drell-Yan process is usually parametrized as 



1 da 
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where 9 and <j) are, respectively, the polar angle and the az- 
imuthal angle of dileptons in the Collins-Soper frame [22] . 
At leading order, the (^-independent term of the unpolar- 
ized Drell-Yan cross section for proton-antiproton colli- 
sions is 



da*? 



r (l + cos 2 9) 



df2dx 1 dx 2 d 2 q T 12 M 2 

x ^2 e 2 I d 2 k T d 2 p T 5 2 (k T + p T - q T ) 

a •* 



(2) 



where M 2 is the invariant mass squared of the lepton pair 
and qT is the transverse momentum of the virtual photon 
in the frame where the two colliding hadrons are collinear. 
In Eq. © we denoted by fi(x, k T ) the unintegrated quark 
number density in the proton, and we used the fact that 
the antiquark distributions in the antiproton are equal to 
the quark distributions in the proton. 

We will consider the v term in |T]). It is known that 
gluon radiation processes give rise to a non-zero cos 2<j) 
asymmetry, which in case of qq annihilation dominance 
is given by v = Q T /(M 2 + 3Q T /2) [23]. Various quan- 
titative analyses [24,25,4] show that perturbative correc- 
tions are unable to reproduce both the magnitude and the 
Qr-dependence of v as observed by NA10 in the region 
M ~ 4 — 8 GeV, but for lower M values the perturbative 
asymmetry might be relevant. In the present paper we 
will focus on the Boer-Mulders contribution to the cos 2<f> 
asymmetry, which has been shown to explain the NA10 re- 
sults [2] and represents a large effect in the moderate-Qr 
region we are interested inQ 



1 Concerning the perturbative contributions to azimuthal 
asymmetries, a forthcoming analysis of cos 2<j> distributions in 
semiinclusive DIS shows that at moderate Q 2 and low Pt the 
Boer-Mulders effect largely dominates over the perturbative 
contribution (at large Q 2 the situation is reversed) [26] . 



The expression for the contribution of the Boer-Mulders 
distribution to the unpolarized cross-section is [5j 
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where h = qr/Qr, with Qt = |qr|- This term, with its 
peculiar angular dependence, gives a parton model expla- 
nation for the cos 2<f> asymmetry observed in the unpolar- 
ized irN Drell-Yan process [4]. 

From Eqs. ^ and §5§ we get the following expression 
of the coefficient v in (fTJ) (with A = 1, fj, = 0) : 



_ ^ a e 2 a n[hi\hj a ] 

E«eS^[/i ,/i ] ' 
where we used the notations: 

•^A a i /i] = ^2 / d " k T d 2 p T S 2 (k T + p T q T ) 
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H{hi\ hi a ] = ~^2 J d 2 k T d 2 p T S 2 (k T + p T - q T ) 
(2h • krh • Pt — kr • Pt) 
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xhi a (x 1 ,k 2 T)hi a (x2,p 2 T ) 
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2tt 



k T + Qt ^ cos x — 2 k T cos 2 x 

x 2 
m N 

xhi a ( Xl ,k T )hi a (x2,\<lT-kT\ 2 ), (6) 

and x is the angle between qT and kT. The asymmetry 
coefficient v given in Eq. ((4| depends on the kinematic 
variables xi, X2, M and Qt- The GSI-HESR kinematics 
probes the large x region, where the valence quarks and 
in particular the u sector dominate, so that Eq. ((4J can 
be simplified to: 



2H[hi u ,hj u } 



(7) 



Another promising process to study transverse spin 
physics at GSI-HESR is J/ip production, which has been 
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proposed [27] as a method to access transversity by mea- 
suring double spin asymmetry. The J '/tp production events 
are identified by a peak in the dilepton invariant mass 
spectrum at M — mji^ = 3.097 GeV, so one can choose 
M 2 ~ 9 GeV 2 . At the GSI-HESR energy scale (s = 30-45 
GeV 2 in the fixed target mode, or s = 100 - 200 GeV 2 
in the collider mode) J/ip production is dominated by qq 
annihilation [5H]. Since the J/ip is a vector particle and 
the qqJ/ip coupling has the same helicity structure as the 
qqj* coupling, one can get the J/ip production cross sec- 
tion by replacing the quark electric charges by the J/ip 
vector couplings to qq and to , 



1^22 2 

16tt a cm e a 



(8) 



and the virtual photon propagator by a Breit-Wigner func- 
tion, 

J_ _ I ( 9 ) 

where Fj is the J/ip width. This model successfully ac- 
counts for the SPS J/ip production data at moderate s 

m ■ 

Therefore the unpolarized cos2</> asymmetry in the 
J/ip resonance region reads 



(10) 



with only e 2 replaced by (gX) 2 m Eq.(|3]). Since the u quark 
dominates, this asymmetry reduces to Eq. ([7]), which is 
the same as for continuum Drell-Yan, with M 2 ~ 9 GeV 2 . 
Integrating over x\ and x 2 we get 



u(Q T ,M) 



2jd Xl J dx 2 H[hi u , /if 11 ] (5(1 - r 
fdx 1 fdx 2 TlKJ?}S(l-T) 



(11) 



where r = M 2 / x\x 2 s. Notice that the relation 1 — r = 
is valid for Q\ <C M 2 . The rapidity is defined as y — 
g ln(^r) and Feynman's variable is Xp — X\ — x 2 . 



3 Calculation of the cos 20 asymmetry 

To evaluate the cos 2cp asymmetry given in (0) one needs 
to know the form of the fey-dependent distributions ap- 
pearing in the transverse momentum convolution. Useful 
information on the Boer-Mulders function /if (x, ky) can 
be obtained from the study of the cos 2<p azimuthal asym- 
metry in the unpolarized ttA^ Drell-Yan processes, which 
has been measured by the NA10 Collaboration [4] and 
the E615 Collaboration [5]. In pUlfTT] this asymmetry was 
estimated by computing the /if distribution of the pion 
and of the nucleon in a quark spectator model [T2f30] and 
was compared with NA10 data. To compute the cos 20 
azimuthal asymmetry in the unpolarized pp collision we 
adopt the same distributions /if (x, k 2 -,) and fi(x, k T ) used 
in [TT] . We assume that the observables are dominated by 
u quarks. The set of the transverse-momentum dependent 



distribution functions is (for simplicity, we consider a spec- 
tator scalar diquark [TT|ll2j ) 

= *(!-*)» (l 7 L tTkj.V +14 ' < 12 > 
t> , 4 ,4, |1 -,,. ;<- : ^ , ( ,3, 

where N is a normalization constant, m q is the constituent 
quark mass, and 

L 2 = (1 - x)A 2 + xm 2 d - x (1 - x)m 2 N . (14) 

Here yl is a cutoff appearing in the nucleon-quark-diquark 
vertex and md is the mass of the scalar diquark. As it is 
typical of all model calculations of quark distribution func- 
tions, we expect that Eqs. (fT2|) and (fT3|) should be valid at 
low Q 2 values, of order of 1 GeV 2 . The average trans- 
verse momentum of quarks inside the nucleons, as com- 
puted from [fl2"]). turns out to be (fc^) 1/2 ~ 0.54 GeV. For 
the parameters in Eqs. (fT2"|) and [T3"l) , we choose the values 
m d = 0.8 GeV, m q = 0.3 GeV, A = 0.6 GeV, a s = 0.3, 
which are the same as in [TT | I14 | . 

Compared with previous calculations of Drell-Yan asym- 
metries [9, 18 , our model differs for the explicit form of the 
proton-quark-diquark effective coupling and for the values 
of the parameters. We adopt a dipole form factor [T^] for 
the effective coupling, whereas in other computations this 
coupling is taken as a constant [9], or as a Gaussian form 
factor [18] . This leads to different predictions on the mag- 
nitude and the (Jr-dependence of the cos 2cp asymmetry, 
as we will see below. Our model proved to be capable to re- 
produce the observed cos 2<p asymmetry of the NA10 7r~ 
Drell-Yan data with the correct Qr-dependence [11] . 

The Drell-Yan and J/ip production events are recon- 
structed from the dilepton invariant mass spectrum, where 
the J/ip events correspond to the peak at the invariant 
mass M ~ Mj/^ ~ 3 GeV, and the genuine Drell-Yan 
events correspond to the continuum spectrum below and 
above the J/ip peak. For the experimental programs at 
GSI-HESR, the center of mass energy is s = 30 - 45 GeV 2 
in the fixed target mode and s — 100 — 200 GeV 2 for the 
collider option. With s = 45 GeV 2 and s = 200 GeV 2 
one has on the J/ip peak t = xix 2 — M 2 /s ~ 0.2 and 
t ~ 0.05, respectively. In these kinematic domains va- 
lence quarks dominate and the assumptions made above 
are justified. 

The cos 2(p asymmetry v as a function of Qt is shown 
in Fig. Q] for the fixed-target GSI kinematics. For a quan- 
titative comparison, in the same figure we present the pre- 
diction of a model with constant quark-spectator coupling 
[9], obtained using the same parameters of [9] but for the 
kinematic domain we are considering here. We notice that 
both the magnitude and the Qt dependence arc quite dif- 
ferent, similarly to what happens in the 7rA^ Drell-Yan 
process [11] . The shape of v in the model with a Gaussian 
quark-spectator coupling [TS] also differs from ours and is 
similar to that of [9]. 

The Q-r-dependence of the cos 2<p asymmetry v in J/ip 
production is shown in Fig. [5J As one can see, the asym- 
metry is sizable and increases with Qt- Some difference 
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Q T GeV 

Fig. 1. The cos 2<p azimuthal asymmetry of Drell-Yan produc- 
tion in pp collisions as a function of Qt, for s = 45 GeV 2 and 
M 2 = 5 GeV 2 . The solid curve is our prediction. The dashed 
curve is the result of [9] extended to the kinematic region of 
interest here. 
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Fig. 2. The cos 20 azimuthal asymmetry of pp —> J/ip X — > 
l + l~X as a function of Qt, for two values of s in the GSI-HESR 
kinematic domain. 



between the two curves (corresponding to two values of s) 
emerges at large Qt, where the asymmetry is suppressed 
for the larger center of mass energy. In Fig. [3] we plot v 
as a function of xf = X\ — xi- The asymmetry lies in the 
range 0.2-0.3. 

J I ip production has been recently suggested by Ansel- 
mino et al. |27] as a candidate reaction to probe the nu- 
clcon transversity by measuring the double transverse spin 
asymmetry in polarized p^ffi process. The dilepton pro- 
duction rate on the J/ip peak is two orders of magnitude 
larger than in the continuum region above the J/ip mass. 
This makes J/ip at GSI energies an excellent process to 
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Fig. 3. The cos 20 azimuthal asymmetry of pp — > J/ipX — > 
l + V X as a function of xf, for two values of s in the GSI-HESR 
kinematic domain. 



measure transverse polarization asymmetries. Here we ap- 
plied the same idea to the unpolarized case, and in partic- 
ular to the cos 2(f) asymmetry. Given the large value of v 
that we found, we can say that the experimental study of 
J/ip and Drell-Yan production from unpolarized proton- 
antiproton collisions at GSI represents a very important 
source of information on the Boer-Mulders distribution 
(of course, a complete analysis of these processes will also 
require a careful consideration of the perturbative effects) . 

4 Conclusion 

In summary, we calculated the cos 2(p azimuthal asym- 
metry of the unpolarized pp dilepton production in the 
continuum region and on the J/ip peak for the GSI kine- 
matics, relying on a model of the Boer-Mulders function 
previously adjusted to fit the available experimental data 
on ttN Drell-Yan. The asymmetry turns out to be rather 
large, of order 0.2-0.3. The size of the asymmetry and 
the high counting rate in the J/ip resonance region make 
the dilepton production in pp scattering at moderate ener- 
gies a very promising process to detect the Boer-Mulders 
function in the early unpolarized stages of the future GSI 
experiments. 
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